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Abstract: We present results of time evolution of oblique MHD plane shocks including diffusive cosmic 
ray acceleration with backreaction on the plasma flows. The simulations include self-consistent effects 
of finite Alfven wave propagation and dissipation. From the computed cosmic ray particle phase space 
distributions we calculate expected leptonic and hadronic emissions resulting from interactions between 
the cosmic rays, magnetic fields, the thermal particle population and relevant astrophysical photon fields. 
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Introduction 

Cosmic ray acceleration in strong shocks is highly 
efficient and naturally leads to substantial modi- 
fication in the shock structure compared to ordi- 
nary gas or MHD shocks described by Rankine- 
Hugoniot relations. The modified shock compres- 
sion can greatly exceed that of an adiabatic gas 
shock and the shock structure includes a foot or 
precursor where plasma is compressed and heated 
as it flows against cosmic rays (CRs) streaming 
ahead of the relatively thin dissipative shock tran- 
sition. Since the cosmic rays CRs are accelerated 
by diffusive propagation through the entire shock 
structure, these shock modifications also alter the 
spectrum of the CRs compared to the test particle 
spectrum formed in a discontinuous transition. 

The presence of magnetic fields is essential to the 
physics of diffusive shock acceleration (DSA), be- 
cause the principal CR scattering mechanism is gy- 
roresonant interaction with Alfven waves. That is 
typically modeled by way of the spatial diffusion 
coefficient for the CRs in an otherwise gasdynamic 
model of the shock (e.g., [3, 4]). In some calcula- 
tions the influence of finite streaming of the Alfven 
waves with respect to the bulk plasma and the lo- 
cal dissipation of wave energy (i.e., "Alfven trans- 
port") have been included (e.g., [1, 7, 5]). On the 
other hand, despite the facts that typically one ex- 
pects an oblique magnetic field with respect to the 
shock normal and that the inclusion of significant 



Alfven transport effects imply significant MHD ef- 
fects, very few calculations have been published 
that include a full MHD treatment of DSA, espe- 
cially when the CR spectrum is calculated self- 
consistently [8, 11]. To explore the importance of 
such a self-consistent treatment we include all of 
the above physics in the calculations reported here. 
Since modified CR shocks are generally evolving 
structures so long as the CR spectrum continues 
to extend to higher energies, our treatment is also 
time dependent. 

The CRs accelerated in astrophysical shocks 
will produce observable electromagnetic emissions 
through their interactions with the local plasma and 
ambient radiation fields. The intensity and spec- 
tra of these emissions will generally depend on 
the CR spectrum as well as the structure of the 
shock. To illustrate the importance of MHD ef- 
fects in these emissions we include calculations of 
emissions generated by both leptonic and hadronic 
CRs in the modified MHD shocks we present here. 

Methods and Model Parameters 

We carried out our simulations in one spatial di- 
mension, x, using our "Coarse Grained finite Mo- 
mentum Volume" (CGMV) scheme for solving 
the CR diffusion-convection equation [9]. The 
CGMV scheme evolves the first two momen- 
tum moments of the CR momentum distribution 
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Figure 1: Top: Immediate-postshock CR proton 
spectra at the gas subshock in four MHD shocks 
as outlined in the text and in Table 1. Bottom: 
Intensity of EM emissions from the computed 
shocks including synchrotron, inverse Compton, 
bremsstrahlung and secondary pion decays. Num- 
bers correspond to models in Table 1 . 

function, f(t,p,x), over finite momentum bins, 
Amp ~ 1, assuming a piecewise powerlaw mo- 
mentum dependence for f(t,p,x). The power- 
law slope in each bin is part of the obtained so- 
lution. We have demonstrated that this approach 
provides accurate solutions to the dynamics and the 
evolution of the CR distribution at greatly reduced 
computational effort in comparison to finite differ- 
ence methods. The CGMV routines were incorpo- 
rated into our well-tested TVD MHD code [6] and 
CR modified shock solutions were tested against a 
chain of previously published simulations. 



Model 


M AyX 


#0 


S (mG) 


x D , (cm) 


1 


3000 





0.46 


2.3(14) 


2 


30 





46.0 


2.3(12) 


3 


30 


30 


53.0 


2.0(12) 


4 


30 


75 


180.0 


5.7(11) 



Table 1 : Shock Model Parameters 



We assume upstream of the dissipative subshock 
that the net scattering turbulence velocity with re- 
spect to the plasma is the Alfven velocity paral- 
lel to the local magnetic field vector in response to 
resonant amplification of upstream-facing Alfven 
waves by streaming CRs. Downstream of the sub- 
shock we assume isotropic Alfvenic turbulence. 
We assume also that the Alfvenic turbulence dis- 
sipation rate matches the local growth rate pro- 
duced by resonant scattering; i.e., — va ■ VP C = 
—VAxdPc/dx [1]. The simulations include both 
a dynamically important CR proton and a passive 
CR electron component. Both components have 
the same spatial diffusion model; the electrons dif- 
fer in their evolution only through their energy 
losses to synchrotron emission and inverse Comp- 
ton emission of a combined cosmic microwave 
background and galactic interstellar radiation field. 
Spatial CR diffusion is isotropic and "Bohm-like 
", resulting in an effective diffusion coefficient 
along the shock normal, k = kq(Bq/ B)p, with 
kq = (l/3)(m p c 3 /(eE>o). We henceforth express 
all particle momenta in units of m p c. Subscripts 
refer to conditions far upstream. 

Initially the upstream CR population is void. CR 
protons and electrons are injected at the subshock 
following the simple presription that a fraction ei n j 
of the thermal particle population passing through 
the shock is injected into the CR population with 
momenta pi n j — m p Xc s .2, where c s .2 is the post- 
shock sound speed, and we set A = 2. In the sim- 
ulations presented here the proton injection frac- 
tion is €i n j :P — 10 -2 , while the electron value is 

f- — io'~ 4 
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These shocks all have initial sonic Mach num- 
bers, M Sj o = u s /csfl — 30, with physical shock 
speed, u s = 0.01c = 3000 km/s. The large 
scale magnetic field is placed in the x — y plane. 
We include two parallel shocks (B y = 0, 9 = 
tarT^Byfil B x = 0) with Alfvenic Mach num- 
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bers, Ma = u s /va x = 3000, and 30. There 
are also two oblique MHD shocks with Max = 
u s /vax = 30 and 6 = 30°, and 75°. The up- 
stream plasma density, po = Tn p riQ = 1-67 x 
10~ 24 g/cm 3 . For a given Alfvenic Mach number 
this fixes the upstream magnetic field strength, as 
given in Table 1 . Models 2-4 nominally all have the 
same Alfven wave advection and dissipation prop- 
erties. 

Results 

The properties of the four simulations at t = 180 
are presented in Figures 1 and 2. By this time 
each of the shocks is close to reaching its asymp- 
totic postshock compression and pressure values, 
although the shock transition and the CR spectra 
would continue to spread self-similarly [9]. The 
characteristic diffusion length, xd,o = kq/u s , and 
diffusion time, to,Q — ^o/it 2 , along with the up- 
stream mass density, po provide convenient scal- 
ing units. Setting ko = 1, we have u s = 1 
and an upstream gas pressure, P g ,o — 1/1500. 
The magnetic field is presented in units such that 
the magnetic pressure, Pb = (1/2)B 2 . The full 
computational box in each case spans the spatial 
domain [—100, 100], with the associated physical 
x D fl listed in Table 1. 

The four shocks are shown at equivalent times in 
the sense that test particle DSA theory would pre- 
dict the same maximum CR momentum, p max ~ 
(l/8)t ~ 20 (E max ~ 20GeV) (e.g., [10]). While 
that prediction is roughly confirmed in Figure 1, 
it is also clear that the particle distributions, the 
predicted emissions and the shock structures show 
obvious differences among the models. To under- 
stand the differences it is easiest to begin with a 
comparison of the shock structures as illustrated in 
Figure 2. Here we see that the total shock compres- 
sion is significantly reduced in all the models with 
Max = 30 in comparison to the Ma x — 3000 
case. That is mostly the result of increased Alfven 
wave dissipation in the shock precursor and a re- 
duced net velocity change sensed by CRs across 
the shock (e.g., [10]). That also leads to a reduced 
efficiency in CR acceleration, as pointed out by 
a number of previous authors (e.g., [2]). There 
is a further reduction in shock compression and 



DSA efficiency when the magnetic field is oblique, 
because the transverse magnetic field component 
contributes a significant pressure gradient that re- 
sists compression through the shock. The total 
shock compression in the 6 = 75° model 4 is 
only 1/3 that of the parallel, essentially gasdynamic 
model 1, and the postshock P c is reduced by about 
a factor of two in the same comparison. The par- 
ticle spectra respond to these trends through a re- 
duction in the concavity below p max , since the re- 
duced compression reduces the spread in velocity 
changes sensed by particles as they scatter through 
the precursor. We note for these model parameters 
and simulation times that the electron and proton 
CR distributions do not differ substantially. 

Although the differences in particle spectra exhib- 
ited in Figure 1 are relatively small, they translate 
into substantial differences in the expected electro- 
magnetic spectra. This is illustrated in the right 
panel of Figure 1, where we present the inten- 
sity of radiation from synchrotron, inverse Comp- 
ton, bremsstrahlung and secondary pion decay pro- 
cesses found by integrating along a line of sight 
parallel to the shock normal. The listed order of 
processes corresponds to the order of dominance in 
features seen in the intensity plot beginning at low 
frequencies. The large variation in synchrotron in- 
tensity reflects the two orders of magnitude range 
in magnetic field. More interesting, perhaps, is 
the large range in the bremsstrahlung and pion de- 
cay emissions coming from the large reductions 
in shock compression and DSA efficiency in the 
MHD shocks. 

Conclusions 

MHD CR shocks evolve in noticeably different 
ways in comparison to gasdynamic CR shocks. Fi- 
nite Alfven speeds reduce the efficiency of diffu- 
sive shock acceleration. Magnetic pressure gra- 
dients through the shock transition reduce com- 
pression, further reducing acceleration efficiency. 
These effects can substantially alter predictions 
of the nonthermal emissions associated with the 
shocks, in particular reducing nonthermal X-ray 
and 7-ray emissions. 
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Figure 2: Profiles of the four MHD shocks. Left to right: Gas density, p, and pressure, P g , CR pressure, P c , 
and the transverse magnetic field, B y . Numbers in the P c plot correspond to models in Table 1. 
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